Non-Markovian quantum dynamics: local versus non-local 
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We analyze non-Markovian evolution of open quantum systems. It is shown that any dynamical 
map representing evolution of such a system may be described either by non-local master equation 
with memory kernel or equivalently by equation which is local in time. These two descriptions are 
complementary: if one is simple the other is quite involved, or even singular, and vice versa. The 
price one pays for the local approach is that the corresponding generator keeps the memory about 
the starting point 'to'. This is the very essence of non-Markovianity. Interestingly, this generator 
might be highly singular, nevertheless, the corresponding dynamics is perfectly regular. Remarkably, 
singularities of generator may lead to interesting physical phenomena like revival of coherence or 
sudden death and revival of entanglement. 
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The non-Markovian dynamics of open quantum sys- 
tems attracts nowadays increasing attention It is 
very much connected to the growing interest in control- 
ling quantum systems and applications in modern quan- 
tum technologies such as quantum communication, cryp- 
tography and computation It turns out that the 
popular Markovian approximation which does not take 
into account memory effects is not sufficient for modern 
applications and todays technology calls for truly non- 
Markovian approach. Non-Markovian dynamics was re- 
cently studied in (sl-fisj. Interestingly, several measures 
of non-Markovianity were proposed during last year 

The standard approach to the dynamics of open sys- 
tem uses the Nakajima-Zwanzig projection operator tech- 
nique (20| which shows that under fairly general condi- 
tions, the master equation for the reduced density matrix 
p{t) takes the form of the following non-local equation 



IC(t ~ u)p{u) du , p(to) = Po , (1) 



in which quantum memory effects are taken into account 
through the introduction of the memory kernel IC(t): this 
simply means that the rate of change of the state p{t) at 
time t depends on its history (starting at t = to). Usu- 
ally, one takes to = 0, however, in this letter we shall 
keep '^o' arbitrary. An alternative and technically much 
simpler scheme is the time-convolutionless (TCL) projec- 
tion operator technique [H, Ull, in which one obtains 
a first-order differential equation for the reduced density 
matrix. The advantage of the TCL approach consists 
in the fact that it yields an equation of motion for the 
relevant degrees of freedom which is local in time and 
which is therefore often much easier to deal with than 
the Nakajima-Zwanzig non-local master equation ([T|). 

An essential step to derive TCL from ([T]) relies on the 
existence of certain operator inverse [22!]. However, this 
inverse needs not exist and then the method does not 
work [Ull^]. Moreover, even if it exists the correspond- 
ing local in time TCL generator is usually defined by the 



perturbation series (see e.g. detailed discussion in [l|) in 
powers of the coupling strength characterizing the sys- 
tem. However in general the pcrturbative approach leads 
to significant problems. For example the dynamical map 
needs not be completely positive if one takes only finite 
number of terms from the perturbative expansion. 

In the present paper we take a different path. We 
show that any solution of the non-local equation ([l} al- 
ways satisfies equation which does not involve the integral 
memory kernel, i.e. it is local in time. However, the cor- 
responding generator is effectively non-local due to the 
fact that it keeps the memory about the starting point 
to. Moreover, as we shall see, this generator may be sin- 
gular, nevertheless, it always leads to perfectly regular 
dynamics. 

Let us start with the standard Markovian master equa- 
tion 



dpjt) 
dt 



Cp{t) , p{to) = Po , 



(2) 



where £ is a time- independent generatorpossessing the 
following well known representation j23l . [24 | 



Cp = -i[H,p\ + ^ (v^pV^ - \MV^,p} 



(3) 



The above structure of C guaranties that dynamical map 
A(t, to), defined by p{t) = A{t,to)po, is completely pos- 
itive and trace preserving for t > to. Note that A{t,to) 
itself satisfies Markovian master equation 



|A(t,to) 



£A(t,to) , A(to,io) = ]l 



(4) 



{t-to)C 



and the solution for A{t, to) is given by A{t, to) = e 
which implies that A{t,to) depends only upon the dif- 
ference 't — to' and hence A{t) := A{t,0) defines a 1- 
parameter semigroup satisfying homogeneous composi- 
tion law 



A(ti)A(i2) = A(ti +t2) , 



(5) 



2 



for < 1,^2 > 0. In general the external conditions which 
influence the dynamics of an open system may very in 
time. The natural generalization of the Markovian mas- 
ter equation ^ involves time-dependent generator C{t) 
which has exactly the same representation as in ([3|) with 
time-dependent Hamiltonian H{t) and time-dependent 
Lindblad operators Va{t). Therefore one gets the follow- 
ing master equation for the dynamical map A(t,to) 

^A(t, to) = C{t) A{t, to) , A(to, to) = t, (6) 
at 

with leads to the following solution 

A(i,to) = Tcxp(^^ Cir)dr^ , (7) 

where T stands for the chronological operator. Clearly, 
A(i, to) no longer depends upon H—to' but it still satisfies 
inhomogeneous composition law 

A(t,s)- A(s,to) = A(t,to) , (8) 

for t > s > to- We stress that ^ although time- 
dependent is perfectly Markovian. 

Let us turn to the non-Markovian evolution ([T]). One 
obtains the following equation for the corresponding dy- 
namical map 

-A{t,to)= / dTlC{t-T)A{T,to) , A(to,io) = ]l- (9) 
dt Jt, 

Now comes an essential observation: A(t, to) does depend 
upon the difference ^t — to'' and hence it shares the same 
feature as the Markovian dynamics with time indepen- 
dent generator ([2|) . The proof is very easy. Observe that 
any non-Markovian dynamics in H may be defined as 
a reduced Markovian dynamics on the extended Hilbert 
space v. ® Ha CHa denotes ancilla Hilbert space) . If lu 
denotes a fixed state of the ancilla, then 

Ait,to)p := Tr,[e(*-*«)^(p®w)] , (10) 

where we trace out over ancilla degrees of freedom and L 
denotes the total Markovian generator in H (E)'Ha- Since 
the r.h.s of ([TU|) depends on 't — to' so does the non- 
Markovian dynamical map A{t,to). Hence, the non- 
Markovian dynamics is homogeneous (depends on t — to) 
but of course does not satisfy the composition law ([5|). 
This is the very essence of non-Markovianity and it docs 
provide the evident sign of the memory. 

Suppose now that A(t,to) satisfies non-local equation 
(in]). Taking into account that A is a function of r — t—to, 
let us consider its spectral decomposition 

A{t)p ^ W F,{T)Tr{Gl{T)p) , (11) 

/J 

where F^{t) and G^{t) define the damping basis for 
A(r), that is, Tt:{F^{t)GI{t)) = 5^,^. Clearly, for t = 
one has Ap(0) = 1. Now one defines the formal inverse 

A~\T)p = Y,\-\T)F,{T)Tr{Gl{T)p) , (12) 



such that A(t)A-1(t) = 1 for r > 0. It should be 
stressed that A~^(r) needs not exist (it does exist if and 
only if Ap(t) 7^ 0). Moreover, the existence of A^^(r) 
does not mean that the dynamics is invertible. Note, that 
even if A~^(t) does exist it is in general not completely 
positive and hence can not describe quantum evolution 
backwards in time. Actually, A~^(r) is completely posi- 
tive if and only if A(t) is unitary or anti- unitary. In this 
case |A^i(r)| ~ 1 and A~^(r) — \f_i{T). It is therefore clear 
that the non-local equation (|9]) reduces formally to the 
following one 

^A{t,to)=C{t~to)A{t,to) , A{to,to) = l, (13) 
dt 

where the time-dependent generator C{t) is defined by 
the following logarithmic derivative of the dynamical map 

£(r):^|-A(r).A-i(r). (14) 
One easily finds the following formula 

/:(r)p = ^£^,(r)Tr(Gt(r)p) , (15) 

with 




where for simplicity we omit the time dependence. In 
particular, if the damping basis is time- independent, and 
A^(t) = e'^'f then the spectral decomposition of C{t) 
has a particulary simple form 

C{T)p = J2%ir)F,Tr{Glp) . (16) 

Summarizing, we have shown that each solution A{t,to) 
to the non-local non-Markovian equation (O docs satisfy 
the first order differential equation ([T^ . Let us observe 
that Eq. is local in time but its generator does re- 
member about the starting point 'to'. This is the most 
important difference with the time-dependent Markovian 
equation The appearance of 'to' in the generator 

£(t — to) implies that C is effectively non-local in time, 
that is, it contains a memory. Therefore, the local equa- 
tion p3p is non-Markovian contrary to the local equation 
^ which does does not keep any memory about to. Note, 
that solution to (fT^ is given by 

A(t,to) = T exp £(r)dT^ . (17) 

It shows that A(t, to) is indeed homogeneous in time (de- 
pends on 't — to'). However, contrary to (O, it does not 
satisfy the composition law. Again, this is a clear sign 
for the memory effect. 

One may ask a natural question: how to construct 
non-Markovian generator C{t). The general answer is 
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not known but one may easily propose special construc- 
tions. Let £ be a Markovian generator defined by Q and 
define C{t) — a{T)C. It is clear that if /J" a{u)du > 
for r > 0, then A(r) = exp{J^ a{u)du C) defines com- 
pletely positive non-Markovian dynamics. This construc- 
tion may be generalized as follows: consider N mutu- 
ally commuting Markovian generators Ci, . . . , Cn and N 
real functions ^^(t) satisfying ak(u)du > 0. Then 
C{t) = ai{T)Ci + . . . + a]^{T)CN serves as a generator of 
non-Markovian evolution. Finally, let us observe that if 
C{t) is a time-dependent Markovian generator (i.e. it has 
the Lindblad form ([3]) with time-dependent Hamiltonian 
H{t) and noise operators Va{t)), then C{t — t^) gener- 
ates the non-Markovian dynamics for t > tQ. We stress 
that these constructions provide only restricted classes 
of examples of non-Markovian generators. All of them 
start with a set of Markovian generators and produce a 
non-Markovian one. It turns out (see Example [3] below) 
that one may construct generators which do not fit these 
classes. 

Let us illustrate our analysis with the following simple 
examples. 

Example 1 Consider the dynamical map for a qudit {d- 
level quantum system) given by 



where 



A(r) 



1 



f{u)du 1+ / f{u)duV , (18) 
J Jo 



where V : B{C^) — > B{C^) denotes completely positive 
trace preserving projection. For example take a fixed 
qudit state w and define V by the following formula Vp = 
wTrp. The real function '/' satisfies: 



< 



f{u)du < 1 , 



for any r > 0. Note that f{u) needs not be positive. 
If f{u) > (u > 0), then A(r) defines quantum semi- 
Markov process and the function f{u) may be interpreted 
as a waiting time distribution for this process 0, [l3| . 
Clearly, A(t) being a convex combination of II and V is 
completely positive trace preserving map and hence it 
defines legal quantum dynamics of a qudit. The corre- 
sponding memory kernel is well known 0, and it is 



given 



/C(t) = k{t)Cq , 



(19) 



where the function k{t) is defined in terms of its Laplace 
transform as follows 



k{s) 



(20) 



and £o is defined by £o = ^ ^ 1- Note, that Cq has ex- 
actly the structure of the Markovian generator ([3]) with 
H — 0, and the Lindblad operators Va define Kraus rep- 
resentation of V, that is Vp = J^a^apV^. One easily 
finds for the corresponding generator 



a(T) 



fir) 



(22) 



Let us observe that 



a{u)du = - In (^1 - f fiu)duj > 



and hence this example gives rise to C(t) = a(r)i2oi with 
Markovian Cq and a(r) satisfying JJ" a{u)du > 0. We 
stress that a{T) needs not be positive. It is positive if and 
only if /(t) corresponds to the waiting time distribution 
[1, [l^. Note the striking similarity between formulae 
(PO)) and (|^^ . It should be stressed that in this case one 
knows an explicit formula for time- local generator C{t). 
Note, however, that in general one is not able to invert 
the Laplace transform of k(s) from the formula ([20|) and 
hence the explicit formula for the memory kernel K,{t) is 
not known. 

Example 2 The previous example may be easily gener- 
alized to bipartite systems. Consider for example a 2- 
qubit system and let T' be a projector onto the diagonal 
part with respect to the product basis |m (g) n) in C^. 
Let us take as an initial density matrix so called X-state 
[2^ represented by 



Po = 



pii pi4 

P22 P23 

P32 P33 

P41 P44 



(23) 



C{t) = a(T)£o , 



(21) 



It is easy to see that A(r) defined by ([THl) does pre- 
serve the structure of X-state, that is, p(t) has exactly 
the same form as in (|23p with r-dependent Pmn- It 
is clear that the diagonal elements are time indepen- 
dent pkkir) = Pkk, and pfe/(r) = (1 - f{u)du)pki, 
for A; 7^ I. The entanglement of the 2-qubit X-state 
p{t) is uniquely determined by the concurrence C(t) = 
2 max{ci(T), C2(t), 0} , where 

ci(t) = |p23(t)|- VP11P44 , C2(t) = |pi4(t)|- VP22P33 , 

that is, p{t) is entangled if and only if ci(t) > 
or C2(t) > 0. Let us observe that the function /(r) 
controls the evolution of quantum entanglement. Con- 
sider for example /(r) = £76"'^'^, with 7 > and 
e e (0,1]. One finds from the following formula 
Q!(r) = £7[(1— e)e''''^-|-e]~"'^. Note, that for £ = 1 it reduces 
to a{T) ~ 7, that is, it corresponds to the purely Marko- 
vian case. Hence, the parameter '1 — e' measures the 
non-Markovianity of the dynamics. Suppose now that 
Pq is entangled. The entanglement of the asymptotic 
state is governed by C(oo) = 2max{ci(oo),C2(oo),0} , 
with ci(oo) = (1 - £)|p23| - VPiipil and 02(00) = 
(1 — £)|pi4| — y/P22P33- It is clear that in the Marko- 
vian case (e = 1) the asymptotic state is always separa- 
ble (C(oo) = 0). However, for sufficiently small '£' (i.e. 
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sufficiently big non-Markovianity parameter '1 — e') one 
may have Ci(cx)) > or 02(00) > 0, that is, the asymp- 
totic state might be entangled. This example proves the 
crucial difference between Markovian and non-Markovian 
dynamics of composed systems. In particular controlling 
'e' we may avoid sudden death of entanglement [25l |. 

Example 3 Consider the pure dccohcrcncc model de- 
fined by the following Hamiltonian H = Hn + Hs + 
HsR, where is the reservoir Hamiltonian, Hs = 
X]„ en-Fri {Pn = |"-)(?^|) the system Hamiltonian and 

i^Si?-^Pn®S„ (24) 

n 

the interaction part, _B„ = being reservoirs operators. 
The initial product state p (g) lor evolves according to the 
unitary evolution e~*^*(p (g) Wfl)e*^* and by partial trac- 
ing with respect to the reservoir degrees of freedom one 
finds for the evolved system density matrix 

p{t) = TrR[e-'"*{p<»ujR)e'"*] = ^ c™„(t)P™pP„ , 

n.m 

where c™„(t) = Tr(e-'^'"*Wfle'^"*), Z„ = e,jR + Hr + 
Bn being reservoir operators. Note that the matrix 
Cmn(i) is semi-positive definite and hence 

A(t) P = Y^ C„,niT)PrapPn ■ (25) 

defines the Kraus representation of the completely pos- 
itive map A(r). The solution of the pure decoherence 
model can therefore be found without explicitly writ- 
ing down the underlying master equation. Our method, 
however, enables one to find the corresponding generator 
C{t). It is given by the following formula 

'C(t) p = ^ a,nn{T)PmPPn , (26) 

where the functions amn(T) are defined by amn = 
Cmn/cmn- It shows that the pure decoherence model may 
be defined by local in time master equation (|13p with the 
non-Markovian generator (pS)) . It should be stressed that 
this generator is not of the Lindblad form. 

Example 4 Consider the non-Markovian dynamics of a 
qubit generated by the following singular generator 

£{T)^t&nT£o, (27) 



with Co being the pure dephasing generator defined by 
CqP = cTzpcTz — P- This generator was analyzed in 
in the context of quantifying non-Markovianity of quan- 
tum dynamics. Note that C{t) has an infinite number of 
singular points t„ = (n -I- i)7r. One easily finds the fol- 
lowing perfectly regular solution for the dynamical map 
A(t) = \{1 + cosT)t + i(l - cobt){Co + 11), that is, the 
density matrix evolves as follows 

p{r)^( Pi2COsr\ 

' \ P21COST P22 J ^ ' 

and hence it displays oscillations of the qubit coherence 
Pi2(t). Note that /o(r„) is perfectly decohered, whereas 
for T = nn the coherence is perfectly restored. Finally, 
one finds extremely simple formula for the corresponding 
memory kernel IC{t) = ■^Cq, for t > to- Hence, one 
obtains (|28p either from the non-local equation with time- 
independent memory kernel IC{t) = 5 'Co, or from time- 
local equation with highly singular generator (j27p . 



In conclusion, we have shown that non-Markovian quan- 
tum evolution may be described either by the non-local 
equation ([ij or by a time- local equation ([T5|) . A sim- 
ilar strategy based on pseudo-inverse maps have been 
recently applied in [26j . We stress, however, that our 
approach is different. Clearly, the local approach is 
technically much simpler, however, the prize we pay for 
this simplification is that the corresponding generator 
C{t — to) is no longer local in time but it contains a mem- 
ory about the starting point Hq'. Our examples show that 
these two descriptions are complementary: if /C(t) is sim- 
ple (like JC{t) = ^ £0), then £(t) is highly singular (like in 
(1^71) ). Vice- versa in the Markovian case C{t) — Cm but 
the memory kernel K. is highly singular and it does involve 
the Dirac delta-distribution /C(r) = 2S{t)Cm- Remark- 
ably, singularities of C might provide interesting physical 
content. Note, that the singularities of 'tanr' in Exam- 
ple |4] imply the interesting features of the dynamical map 
if we evolve a maximally entangled state P+ of two 
qubits via the channel ^{t) := 1(8) A(t), then 4'(t)P+ 
is separable if and only if r = r„. It shows that the dy- 
namics 4'(t) gives rise to entanglement sudden death [23 
whenever C(t) is singular and then entanglement starts 
to revive. 
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